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We consider quantum corrections to classical real time correlation functions 
at finite temperature. We derive a semi-classical expansion in powers of Ti 

S3 



Oh 



with coefficients including all orders in the coupling constant. We give ex- 
plicit expressions up to order h 2 . We restrict ourselves to a scalar theory. 
This method, if extended to gauge theories, might be used to compute quan- 
tum corrections to the high temperature baryon number violation rate in the 



X 

5— i ' Standard Model. 



I. INTRODUCTION 

In the Standard Model baryon number is not conserved JT|. While the rate for baryon 
number violating processes is exponentially suppressed at zero temperature it can be signif- 
icant at temperatures of the order of the electroweak scale This could have had an 
influence on the baryon asymmetry of the universe (for a recent review see || ) . 
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Baryon number nonconserving processes are mediated by topologically nontrivial 
vacuum-to- vacuum transitions in the electroweak theory. These transitions are characterized 
by a change AiVcs = ±1 of the the Chern-Simons number 

Ncs = ^e ijk J £x (f$AI - 9 -e ahc A«A)Al} (1) 

where A® are the SU(2) gauge fields and g is the weak coupling constant. The change of 
baryon number AB is related to AiVcs by 

AB = n f AN cs (2) 

where rif is the number of fermion families. At finite temperature the iVcs _ changing transi- 
tions occur by thermal fluctuations across an energy barrier. For temperatures less than the 
critical temperature of the electroweak phase transition T c ~ 100 GeV the transition rate 
can be calculated in a weak coupling coupling expansion around the the so called sphaleron 
which is a saddle point on this barrier. 

The sphaleron approximation can not be extended to temperatures above T c . The field 
configurations which dominate the rate for T T c have a typical size of order ~ 1/ {g 2 T) p] . 
This is the scale at which perturbation theory at finite temperature breaks down ||. For this 
temperature range the baryon number violation rate per unit volume has been estimated as 

1 

T = K{a w Tf (3) 

where k is a dimensionless constant and aw = g 2 /(47r). 

Grigoriev and Rubakov [|10j have suggested a numerical algorithm for computing T. 
They argued that the thermal activation across a barrier is a classically allowed process 
and that one should determine the rate in a purely classical field theory. This consists of 
solving the classical equations of motion for the gauge fields and measuring the change of 
the Chern-Simons number. Finally the average over initial conditions has to be performed 
with the Boltzmann weight exp(—(3H). Since the relevant field modes are long wavelength 
and involve many quanta one would expect that they in fact behave classically. 
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The quantity of interest is the real time correlation function 

Ccs(t) = (Wcsit) - iV cs (0)] 2 ) (4) 



where (■ ■ •) denotes the average over a thermal ensemble. It is commonly assumed that the 
Chern-Simons number performs a random walk so that for large times Ccs{t) grows linearly 
with t. Then the rate T is related to Ccs{t) through 

Ccs(*) - vrt (5) 
for t —>■ oo where V is the space volume [0. Using real time evolution of classical equations of 



motion on a Hamiltonian lattice Ambj0rn et al. |TT] obtained a lower bound for T. Ambj0rn 
and Krasnitz [JTJ] performed a more precise calculation for the pure SU(2) gauge theory, they 
found T = l.l(awT) 4 . This result has been confirmed by Moore [|13| , |14"f who considered the 



change of Ncs in the presence of a chemical potential for Nqs- Tang and Smit |15| have used 
the algorithm of Ref. |T^] to study the the SU(2)-Higgs model. For T ^> T c their results 
agree with those of Ref. [E|. For temperatures less than T c they found a result which is 



much larger than the one obtained in the sphaleron approximation 0,0. Moore and Turok 
[ rS| found similar results for T < T c but they argue that for this temperature range the 



classical approxomation for the correlation function might not be reliable. 

There are several questions arising in this context. It has been known since the work of 
Rayleigh, Jeans, Planck and Einstein that statistical mechanics for a field theory is ill defined 
due to ultraviolet divergences. Many physical quantities depend on the cutoff and do not 
have a finite continuum limit when computed classically. For static quantities the classical 



theory is equivalent to a three dimensional Euclidean quantum field theory p7|-pT|. It is 
possible to introduce counterterms in the classical Hamiltonian such that the continuum limit 
exists. Since the three dimensional theory is super-renormalizable only a finite number of 
divergences occur. For time dependent quantities the situation is more complicated. In gauge 
theories there are additional ultraviolet divergences not related to the three dimensional 



theory |2^| . These arise in diagrams which in the corresponding quantum theory give rise to 



terms oc T 2 , the so called hard thermal loops p3| . These are explicitly time dependent. The 



hard thermal loops are due to large loop momenta k ~ T which are much larger than the 
momentum scale ~ g 2 T associated with the iVcs-diffusion. If the high momentum modes 
decouple the continuum limit of certain quantities may exist. The results of Refs. [|I1^-|T6| 
indicate that this is the case for T in Eq. ([5]) when computed for temperatures much larger 
than T c . 

If the continuum limit of the classical correlation function for some particular quantity 
exists, the question remains how accurate the classical approximation is. In Ref. [|16| this 



problem has been addressed perturbatively. It would be desirable to find a way to go beyond 
perturbation theory and to find an algorithm for incorporating quantum corrections into the 
classical lattice calculations done so far. 

In the present paper we address the latter problem. We develop a semi-classical expansion 
of real time correlation functions at finite temperature starting from the quantum mechanical 
correlator. We keep Planck's constant h 7^ 1. We derive an expansion in powers of h. In 
the limit h — > we recover the classical limit. We consider correlation functions of the type 

C(t) = (^[A{0)B{t) + B(t)A(0)]) (6) 

where A(t) and B(t) are functions of the Heisenberg field operators at time t. The operator 
ordering in Eq. (|6]) is the relevant one for the Chern-Simons number diffusion. The problem 
of ultraviolet divergences related to real time evolution will not be addressed here. 

The semiclassical expansion for static finite temperature correlation functions has been 



developed in the early days of quantum mechanics [24] and can be found in textbooks, 



e.g., Ref. The author is not aware of such an expansion for time dependent quantities 
existing in the literature. 

In Section [TJ we illustrate the method for the simplest quantum mechanical system, i.e., 
for one degree of freedom. The generalization to a system with an arbitrary number of 
degrees of freedom (Sect. [IT]] ) and thus to field theory (Sect. |TVp is straightforward. Finally 
we summarize and discuss the results in Sect. |V|. 



II. QUANTUM MECHANICS WITH ONE DEGREE OF FREEDOM 



We consider one quantum mechanical degree of freedom q with conjugate momentum p 
and the Hamiltonian H = p 2 /2 + U(q). The potential U is assumed to be a fourth order 
polynomial in q. The quantum operators are denoted by capital letters. Heisenberg and 
Schrodinger operators are denoted by Q(t), P(t) and Q = Q(0), P = -P(O), respectively We 
consider the two point function 



C(t) = Z-Hr - Q(t)Q(0) + Q(0)Q(t) 



-/3H(P,Q) 



(7) 



where Z = tr {exp(— j3H)}. We will find an expansion 



c(t) = c (t) + nc x {t) + h 2 c 2 {t) + ■ 



(8) 



Following |25| we use (Schrodinger) momentum eigenstates \p) to perform the trace and 
insert 1 = J dq\q)(q\ to obtain 

_ x r dpdq 



with 



I=(q\p}(p\^[Q(t)Q(0) + Q(0)Q(t) 



-/3H(P,Q) 



l?>- 



Factorizing the Boltzmann factor one writes I in terms of x as 



(9) 



(10) 



I = exp{-pH(p,q)}x 



(11) 



Differentiating Eq. (|T^) with respect to (3 one obtains p5| 



dpX = -ihp ~ PU'(q)x + d q x 



+ -% 2 [d 2 X - 2f3U'(q)d qX - PU"(q) X + P 2 U'(q) 2 X 



(12) 



where the primes denote derivatives with respect to q. Eq. (Tl2| ) can be solved by iteration 
when x is expanded in powers of h, 
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X = Xo + hXi + h 2 X2 H (13) 

We have to determine the boundary condition at j3 = 0, 

X° = X ((3 = 0) = (q\p)(p\~m)Q(0) + Q{0)Q{t))\q), (14) 
which will be expanded in powers of H as well: 

X° = Xo + hXi + h 2 X° 2 + ---. (15) 



We write x° as 



x °=( q+ l -hd p )r, (10) 



with 

77 = <<7b)(p|Q(t)| 9 ). (17) 

Differentiating with respect to t we find 

^ = ^<gb>(p|[/r(P,g),Q(t)]|g> (18) 

where we have used the fact that the Hamilton operator is time-independent. From Eq. 
( fL8|) we obtain 

dtV=U~ ify + \tfd 2 q + U(q + ihd p ) - U(q))rj 

= pd qV - U'(q)d pV + % -% (d 2 qV - U"(q)d 2 pV ) + h 2 U'"(q)dl V + ^fU^d^. (19) 

The initial condition for 77 is 

rj{t = 0) = q. (20) 

We expand rj in powers of h 

rj = T] + hrix + h 2 rj 2 H (21) 

and solve Eq. (|T9|) by iteration. For tjq we obtain the differential equation 



dtVo = (pd q - U'(q)d p 



= {H,Vo} (22) 
where {, } denotes the Poisson bracket 

{/, 9} = d p fd q g - d p g.d q f (23) 



The general solution to Eq. (gg) is 770 = f (p c (t) , q c (t)) where p c (t) , q c (t) are the solutions 



to the classical Hamilton equations of motion 

q c (t) = {H,q c (t)}, Pc(t) = {H, Pc (t)} (24) 
with the initial conditions 

q B (0) = q, Pc (0)=p. (25) 

With Eq. (g0|) we obtain 

?/o = q c (t) (26) 
and therefore Xo = licit). From Eq. ( |T2"D we see that xo is independent of f3. Thus 

Xo = qq c {t). (27) 

Then we finally obtain 

Co(t) = Z iJ^e-^ qqc (t) (28) 

with Zq = J ^^■e~ l3H ^ p ' q \ which corresponds to the Grigoriev-Rubakov expression. 
Now we consider the order h correction to Eq. (p8|). From Eq. ( |T2| ) we obtain 

Xi = X°i - ip (/3%o - ^ 2 f/'Xo) • (29) 

In the correlator C(t) the function Xi only appears in the form J dpdqexp(—/3H)xi- One 
can integrate this by parts to find the equivalent expression 



Eq. fllED gives 



J dpdqe'^^xi = J dpdqe- pH ^{ X \ ~ \d p d qXQ }. (30) 



Xi = qVi + ^ d pVo (31) 

where 7/1 is determined by 

d tVl = {H, Vl } + l -(d 2 q - U"(q)d*) Vo (32) 

with the initial condition r]i(t = 0) = 0. In order to solve Eq. (|32|) it is convenient to rewrite 
it as 

d t(m ~ = i H , Vi ~ y d P d qVo}- (33) 



Therefore 
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Combining Eqs. fl3"T| ) and ( |3~3| ) we arrive at 



Vi = o3A?c(*)- (34) 



Xi = 2 d P d iXo- (35) 

If now compare Eqs. (|30|) and (|35D we see that 

Ci(t) = 0, (36) 

i.e. the order ft correction to the classical correlation function fl2~8|) vanishes. 

Next we compute the 0(h 2 ) corrections to C . The /^-dependent part of %2 is determined 
by Eq. ([12]). Using the results for xi (Eqs. ( p9|) and (|35|) ) we obtain the lengthy expression 

X2 = X°2 + \ (Pd'xo ~ P 2 U'(q)d qX o - \p 2 U"(q) X o + \^{U'{q)f X ^ 
+ P (l(3d p d? X o ~ \p 2 U'(q)d p d q Xo) 

P 2 (~P 2 d 2 q Xo + \pU"{q) X o + \^U'{q)d qX o ~ \p{U'(q))\ } . (37) 



It can, however, be simplified by partial integrations: 



24 



d 2 q Xo + U"(q)d 2 pX a-PU"(q)xo 



(38) 



We use Eq. ([16]) to write 



X°2 = QV2 + ^pVi 



(39) 



where rj2 satisfies 



d tV2 = {H, m } + l -(d 2 q - U"(q)%)rn + \u"'{q)dl m 



(40) 



with the initial condition f] 2 (t = 0) = 0. 
To compute r] 2 we rewrite Eq. (HO) as 



d t(v2 + \%%Vo) = {H,m + \d 2 p &M + ^U">(q)d 3 p r, 



(41) 



Thus 



V2 = -^d p d q T] + T] 2 



(42) 



where 7/ 2 satisfies 



d t f} 2 = {H,f} 2 }--U'"(q){q,q c (t)h. 



(43) 



Here we have denoted 



{f,g}o = g, {f,g} n +i = {f,{f,g}n} 



(44) 



in order to write the derivatives d p in terms of Poisson brackets. With the initial condition 



772 (£ = 0) we find 



V2 = -^jyt'U'"(q c (t>)){q c (t'), qc (t)} 3 . 



(45) 
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Here we have used the fact that the Poisson bracket is invariant under canonical transfor- 
mations of the variables p and q and that the time evolution can be viewed as a canonical 
transformation. 

Putting everything together we can write the correlation function (|7|) as 

h 2 /3 



C(t) = Z-i I 

+ o(h 3 ) 



■PH(p,q) 



1 + 



24 



(d 2 + U"(q)d 2 - (3U" 



qq c (t) + k qr] 2 



(46) 



where 772 is given by Eq. 



III. MANY DEGREES OF FREEDOM 



The results obtained in the previous section can be easily generalized to a system with 
an arbitrary number of bosonic degrees of freedom. The coordinates qi and momenta pi will 
be collectively denoted by q and p. The notation is completely analogous to the one in Sect. 
[TI]. The Hamiltonian is H = l/2J2iP 2 + U(q). We consider the correlator 



C(t) = Z- [ tr\^[A(Q)B(Q(t))+B(Q(t))A(Q) 



where A(q) and B(q) are some functions of q. We write C(t) as 



The differential equation determining the (3 dependence of x is 25 



fyx = 1 - fop- 



h 



+ y (d 2 iX - 2/3U i (q)d qiX ~ PUu(q) X + P 2 (U t {q)?X 

where a subscript i on U denotes the derivative d qi . 

The boundary condition for Eq. fl49|) is x(P = 0) = x° with 

X° = \(q\p)(p\[A{Q)B{Q{t)) + B{Q{t))A{Q)\\q). 



(47) 



(48) 



(49) 



(50) 
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This can be written as 



X° = -[A(q) + A(q + ihd p )} V 



(51) 



with 



r) = (q\p){p\B(Q(t))\q). 



(52) 



The time dependence of rj is determined by 



9tV = £ (PidqiV ~ Ui{q)d Pi r]) + l -h[Y^d 2 qi r] - £ U i:j (q)d Pi d p , 



1 % 

+ q %2 E u ijk(q)d Pl d Pj d Pk r] + —h 3 Um(q)d Pi d Pj d Pk d Pl rj. 

ijk ijkl 

Note that that this differential equation is completely analogous to Eq. (|l! 
involve the function B which enters only through the initial condition 



(53) 

). It does not 



rj(t = 0) = B(q). 



(54) 



Again we expand rj = t)q + Tvrji + h 2 T] 2 + • • ■ and solve Eq. 
term in this expansion we obtain 



by iteration. For the first 



770 = B(q c (t)) 



(55) 



where q c (t) is the solution to the classical equations of motion with initial condition q c (0) = q, 
q c {0) = p. The 0(h) correction to Eq. (155) is given through 



and the result for 772 is 



(56) 



where 



1 r* 



* = 9 q .d Pi d qj d Pj r] + rj 2 

8 H 



(57) 



f)2 = J dt'U ijk (q c (t')){q cl (t'), {q cj (t'), {q ck (t'), B(q c (t))}}}. 



(58) 
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Now we can write down the first three terms of the expansion flT5|) for x° : 



X°o = A(q)B(q c (t)) 



(59) 



x°i = ^E^Axo 



(60) 



and 



xI = -\T. {^AAAiXl + Aij(q)d Pi d p .B(q c (t))} + A(q)f} 2 



(61) 



with Aij = d qr d qj A. 



Now that we know the boundary condition at /3 — it is straightforward to solve Eq. 
(fl9|) to determine x- The result is somewhat lengthy but it simplifies after integrating 
J Tli(dpidqi) / (2,7ih) exp(—j3H)x by parts. As in Sect. |I|one finds that the 0{Ti) contribution 
to C(t) vanishes and we obtain 



C(t) = Z- 1 [Y[( e-? H ^ 



2nh 



x { 1 + ^ ( E dl + E UiMWvi ~ P E W 



A{q)B{q c {t)) 



ij 



+ o{h 3 



(62) 



IV. FIELD THEORY 



The field theoretical expressions for the correlation function are completely analogous to 
the ones found in the previous section. We consider a theory with one scalar field p and its 
conjugate momentum 7T. We consider the correlation function 

1 



C(t) = Z-\i - $(0, x A )$(t, x B ) + $(t, x B )$(0, x A ) 



-/9H[n,4] 



(63) 



for the Hamiltonian 
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H[ir,<p] = Jd 3 x^ir 2 (x) + U[ V ] 



where 



U[<p]= j d 6 x\-{V V {x)Y + V{ V {x)) 



Replacing 



5 



x- 



Sir(x) Sip(x) 

the h expansion of (|63"D can be read off from Eqs. fl58|) and 



C(t) = Z- 1 J [dird(p]e 



, f« ( 5 2 

1 + -r- / <rx 
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+ J d 3 x'5(x - x") (-V* + V"{<p{x)) 
Scp(x) 2 J 



5(p(x) 2 
5 2 



57r(x)57r(x') 
<p(x A )tp c (t, x B ) + h 2 (p(x A )fj 2 \ + 0{% 3 ) 



where ip c is the classical field satisfying the equation of motion 



<p c - VVc + V'(tp c ) = 



with the initial conditions 



(64) 



(65) 



(66) 



(67) 



(68) 



<f c (0,x) = ip(x), ip c (0,x) = ir(x). 

Furthermore, f\i in Eq. (|B"7D is given by 

r\i = -^l/* dt ' j d 3 xV'''(<f c (t',x)){<f c {t',x),<f c (t,x B )} 3 

where {, }3 is defined as in Eq. ([0]) and the Poisson brackets are 

5f Sg 5g 5f \ 



{f } g} = J d 3 x (- 



(69) 



(70) 



(71) 



Stt(x) S(p(x) Stt(x) 5<p(x) J 

The right hand side of Eq. fl6"7|) contains ultraviolet divergences. Consider, e.g., the term 
5 2 U[(p]/5(p(x) 2 . The singularities arise because we take two functional derivatives of U with 



13 



respect to (p(x), both of them for the same space point x. If we write 5 2 U[ip]/5(p(x) 2 in 
momentum space and introduce a momentum cutoff A it becomes 

Sp-s^H (72) 

The parameters in the Hamiltonian flBlj) are bare quantities which have to be expressed 
in terms of renormalized ones. In the following we will see how the singularities are canceled 
after renormalization in the limit t — > of C(t). In this limit the problem is of course 
much simpler. Ultraviolet divergences related to real time evolution have been discussed 
in Ref. |22] and some of them turned out to be rather nontrivial even in the classical limit 



i.e. without including the 0(h 2 ) corrections. As it was already mentioned the numerical 
calculations done so far indicate that nevertheless the continuum limit of the classical O(h ) 
result for T (Eq. (|5])) exists. Therefore it might even be meaningful to consider quantum 
corrections leaving aside ultraviolet problems for a moment. 
For definiteness we choose the potential in fl65|) as 



V(*0 = ^ 2 + (73) 

For t = the path integral over ir is Gaussian and 7r can be integrated out. This gives a 
constant factor which is canceled by the corresponding factor in Z and which will be ignored 
in the following. Eq. (pTj) becomes 



C(0) =Z~ X J[dtf\e 



/3U[<f] 
2/32 



1 - 



ULJ_ f d*x(-A 5 + -A 3 (m 2 + - V 2 (x)) 
24 2vr 2 J \5 3 V 2 y v '> 



ip{x A )ip(x B ) 



+ ^S(x A -x B )^ + 0(h 3 ). (74) 

We have to express the bare parameters in terms of the renormalized ones which we 
denote by a subscript r. We write 

U = U r + 8U. (75) 
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We include a field independent part in the Hamiltonian which we denote by U°. The 
counterterm 5U° is written as a power series in h, 5U° = 5U® + tfdU® + • ■ ■■ Then one has 
to choose 



§U = _£ J_ f d U I A S + I A 3 m 2 ] (76) 

to cancel the first two singular terms in Eq. (|7^) . 
For the mass renormalization we write 

m 2 = m 2 — <5m 2 . (77) 

We define the renormalized mass in the imaginary time formalism (see, e.g., |[26|| ): 

m r 2 = m 2 + n(a;„ = 0, /i 2 ) (78) 

where II(a; n , /c 2 ) is the self energy for Matsubara frequency u n = n27rT/h and momentum k 
and n is some renormalization point. If fi is chosen sufficiently large, m 2 can be calculated 
in perturbation theory. The 0(A) contribution to Sm 2 is given by the diagram (a) in Fig.|l|: 

2 n= _oo- y ( 27T ) c;2 + k 2 + m 2 
We can expand 5m 2 in powers of h, 

5m 2 = 5ml + h 2 5m 2 2 + 0{h 4 ). (80) 

The zero Matsubara frequency part of Eq. (|79|) gives 

2 ^ (2tt) j k 2 + m 2 

The nonzero Matsubara frequency part can be expanded in powers of h 2 . With a momentum 
cutoff A the first term in this series gives 



Sm} = _Mi_^. (82) 
24 2vr 2 3 



The zero Matsubara frequency parts of all higher loop diagrams for Sm 2 (e.g. the diagram 
(b) in Fig. |l|) contribute to 5m 2 ,. However if one goes beyond O(h ) there are always at least 
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two propagators in which one particular nonzero Matsubara frequency appears. Therefore 
the h expansion for these diagrams beyond O(h ) starts at 0(h 4 ) and Eq. (§2|) is valid to all 
orders in A. For the coupling constant and wave function renormalization there is no 0(h 2 ) 
contribution even at the one loop level. In general there are finite counterterms of the order 
h which, however, do not have to be specified for the present purpose. 



Replacing U — > U T + 5Uo + h 8U2 in Eq. (|74| ) and expanding in powers of h the 0(h 
divergences in Eq. ([74]) cancel and we obtain 



C(0) = Z- 1 / [dtp] exp {-P(U t [cp] + 6U [cp])} 



x^f(x A )<f(x B ) + ^jS(x A - f s ) j + 0(h 3 ) (83) 
This expression is similar to a dimensionally reduced theory |ET| . There are two differences. 



First, the renormalized parameters in Eq. (B3) are different from the ones in Ref. pi 



This difference is due to finite parts in the counterterm 5C/qM corresponding to different 



renormalization conditions. The choice in Ref. 21 was made such that the Greens functions 



in the 3-dimensional theory and the Greens functions with uo n = in the 4-dimensional 



theories agree up to a certain order in the coupling constant. Furthermore, in Ref. |21] the 
renormalized mass was written in terms of the physical Higgs mass. 



The second difference between (|83|) and the expressions given in Ref. |2T] is the appear 



ance of the 5-function in Eq. (|83|). This term arises because in the formalism described 
in the present paper there is no separation of the fields into zero and nonzero Matsubara 
frequency components. Therefore, the correlation function Q53] ) also contains the diagrams 
in which the external lines have nonzero Matsubara frequency. In momentum space the tree 
level propagator for the u n 7^ fields is D(u n , k 2 ) = + k 2 ). Expanding this expression 

in powers of h gives h 2 / (27rnT) 2 + 0(h 4 ). Summing over all n 7^ we obtain the 5-function- 
term in Eq. (|83]). All higher loop contributions to D(cu n , k 2 ) are 0(h 4 ) or higher. Therefore 
they do not contribute to C(0) at the order considered here. 
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V. SUMMARY AND DISCUSSION 



For real time correlation functions at finite temperature we have obtained an expansion in 
powers of Planck's constant H which, at lowest order, corresponds to the classical correlation 



function in the sense of Grigoriev and Rubakov [1C[. We have restricted ourselves to a scalar 
theory. For the correlation functions of the type (||) we found that the first order correction 
in h vanishes. The first nonvanishing contribution is of the order h 2 . The result has been 
written in terms of the solution to the classical field equations and it's functional derivatives 
with respect to the initial conditions. The expansion outlined in this paper can be easily 
extended to higher orders in h. 

The method outlined in this paper might serve as a starting point for numerical compu- 
tations of quantum corrections to classical correlation functions which can not be calculated 
in perturbation theory. An important example is the Chern-Simons number diffusion rate 
in the electroweak theory which determines the rate for baryon number violating processes 
in the early universe. 

The 0(h 2 ) correction to the classical correlation function contains a term ([70| ) which is 
nonlocal in time. One may ask how this term can be implemented in a numerical evaluation 
of the correlation function. It might be possible to find a set of local equations of motion such 
that their solution reproduces Eq. ([7D|). In Eq. floTD the derivatives with respect to the initial 
conditions can be transformed into local operators at t = through partial integrations of 
the path integral. 

An important task is to extend the method outlined here to gauge theories. Another 
important issue are the ultraviolet divergences of the classical result and of the 0(h 2 ) cor- 



rections. For the pure gauge theory considered in Refs. [0,0] this maybe not so urgent 
since there the classical results seem to have a well behaved continuum limit. Then the 
ultraviolet behavior of the 0(h 2 ) corrections may not be that critical either. However, if the 



Higgs field is included as in Ref. [15] the problem of ultraviolet divergences is crucial. 
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FIGURES 




FIG. 1. Feynman diagrams contributing to mass renormalization. 
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